Precision cosmology as a test for statistics 
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We compute the shift in the epoch of matter-radiation equality due to the possible existence of a 
different statistical (non-extensive) background. The shift is mainly caused by a different neutrino- 
photon temperature ratio. We then consider the prospects to use future large galaxy surveys and 
cosmic microwave background measurements to constrain the degree of non-extensivity of the uni- 
verse. 
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A. Introduction 

The possible need of non-extensive thermostatistics in 
cosmology, gravitation and astrophysics has been dis- 
cussed by several authors, for instance see Refs. [0. In- 
deed, the question of to what extent the universe as a 
whole is an extensive system is well justified. It is known 
that whenever long range interactions are present, stan- 
dard (Boltzmann-Gibbs) statistical sums and integrals 
may diverge (see Ref. [|j for a review). In any case, to 
check the validity of a theory we must locate it in a more 
general framework, a group of test theories, which in gen- 
eral will have a free parameter. Then, one has to isolate 
which of the theories of the group are able to reproduce 
experimental or observational facts. Although this is a 
common practice in cosmological models based on alter- 
native theories of gravity, it is not so usual concerning 
alternative statistics. The problem in these cases is try 
to determine if different statistical theories may also yield 
the same observational results. 

The group of test theories will be, for our purposes, 
that introduced by Tsallis in ||. By use of them we 
expect to check the validity of the Boltzmann-Gibbs sce- 
nario. Let us recall the starting points of nonextensive 
statistic (NES), in Tsallis' approach. The formalism be- 
gins by postulating ||: 

Postulate 1.- The entropy of a system that can be found 
with probability pi in any of W different microstates i is 
given by 



S q = 



w / w \ 



(i) 

with q a real parameter. We have a different statistics for 
every possible q— value. In ([!]) we have used, of course, 
that 



(2) 
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Postulate 2.- An experimental measurement of an ob- 
servable A, whose expectation value in microstate i is a%, 
yields the q- expectation value (generalized expectation 
value (GEV)) 



< A 



w 



q 

Pi a i 



Tvp q A, 



(4) 



In general, 



for the observable A. 

These two statements have the rank of axioms. As 
such, their validity is to be decided exclusively by the 
conclusions to which they lead, and ultimately by com- 
parison with observations. We shall not present any fur- 
ther the formalism in this work, referring the reader to 
Refs. [PHq]. A more detailed presentation can also be 
found in [0]. Instead we shall use only the results we 
need for our calculations. This non-extensive formal ap- 
proach was applied in many different systems including 
Newtonian gravity ^,|) , Levy type diffusion M , the prob- 
lem of solar neutrinos JhJ and others. Concerning the 
universe as a whole, we have previously shown that this 
group of test theories will affect primordial nucleosyn- 
thesis yields and the evolution of the neutron to baryon 
ratio in a computable way. Later, this could be used 
to set bounds upon the free parameter they have and 
thus the upon the degree of non-extensivity present in 
the early universe |q, |ll]jn| . In addition, a generaliza- 
tion of the Planck radiation law was derived and applied 
to the cosmic microwave background data (l^Jl^l • In all 
these works, limits upon the non-extensive (freejparame- 
ter were imposed. The typical order of magnitude of the 
bounds is 10~ 4 , what means that possible Bolztmann- 
Gibbs' deviation may be about one part in ten thousands. 

In this letter we shall consider the matter-radiation 
equality epoch. This moment of cosmic history is im- 
portant for every cosmological model because this tran- 
sition changes the growth rate of density perturbations. 
Whereas perturbations inside the horizon are frozen dur- 
ing a radiation dominated universe, perturbations on all 
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scales c an g row (and structure can form) during the mat- 
ter era |n|. Then, the size of the horizon at the time of 
the matter-radiation equality is an important parame- 
ter of the model (indeed, it is the only one for CDM 
models) and it is to be read from the spectrum of den- 
sity perturbations |H| . In the following years, new satel- 
lites will measure cosmic microwave anisotropies and new 
large galaxy surveys will be avalaible. In particular, the 
Sloan Dig ital Sky Survey (SDSS) will acquire 10 6 red- 
shifts |l6| , |l7| . This will make cosmology to enter in an 
era of high precision and about ten cosmological param- 
eters will be accurately determined [Q. In what follows, 
we expect to show that this set of data can also be used 
as a test of the underlying statistical theory To do so, 
we shall compute the shift in the matter-radiation equal- 
ity that would be caused by non-extensivity; in general, 
by the use of a different statistical background. As this 
changes the shape of the power spectrum of density per- 
turbations, something that can be measured from galaxy 
surveys Jl^Jl^l , a new bound upon the statistical descrip- 
tion will arise. In addition, we shall comment on the 
prospect for the accurate determination of the cosmic 
neutrino background and the effective number of neutrino 
species. This will also allows us to extract conclusions on 
which statistical model may be the right description of 
the universe. 



B. Non-extensive corrections and matter-radiation 
equality 

The photon energy density is dominated by the cos- 
mic microwave background. Its temperature is measured 
extremely well p0[ , 



T 7:0 = (2.728 ±0.004)K 
and this gives, for the standard model, 

p 7 ,o = 4.66 x 1CT 34 g cm" 3 . 



(5) 



(0) 



However, if we are to take into account possible non- 
extensive effects, the energy density will in general be 
given by, 

p q = ^9T A + (40.02 ffb + 34.70 ff/ ) T\q - 1), (7) 

where gbj are degrees of freedom of bosons and fermions, 
g = gi, + 7/8gf and q is a real free parameter measuring 
the amount of non-extensivity in the system . When 
q = 1, all results reduce themselves to the standard ones. 

The present radiation density also has a contribution 
from relic neutrinos that, due to the impossibility of de- 
tecting them directly, must be fixed theoretically. The 
relationship between neutrino and photon temperature 
was computed in Ref. Q for a non-extensive framework 
and is given by, 



T — T 



1/3 



(1 + (^-1)0.013). 



(8) 



As we shall see below, equation (g) is an important re- 
lationship, with a number of observational consequences. 
If N v is the number of neutrinos families (hereafter con- 
sidered as three), then 
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^| 3 4.70(g-i; 

7T I 



(9) 



which, adding to p 1 obtained from (Q) leads to, 



Prad = ttL 4 [1 + 0.681 + 10.33(g - 1)] (10) 
15 ' 

for the total radiation density. The first term in the 
bracket is the usual contribution of photons, the second 
stands for the usual contribution of neutrinos while the 
third sums up all non-extensive corrections. 

The Friedmann equation for a zero curvature 
Friedmann-Robertson- Walker universe is 



8ttG 



(11) 



where a is the cosmological scale factor and p is the en- 
ergy density, which has contributions of non-relativistic 
matter and radiation. If these fluids are non-interacting, 
the conservation equation T^ v v = holds for them sepa- 
rately and one gets, 



/a 

Pm — Pm,0 I 

V a 



fa 

Prad — Pr,0 I 

V a 



(12) 



Here, p m ,o stands for the present matter density (which 
is fixed if we assume spatial flatness) and p ra d,o for the 
analogous radiation density. The redshift in the matter- 
radiation equality is given by, 



a 



Pm,0 



Then, 



1 + z, 



eg 



8ttG p 



eq Prad,Q 



24000/j 2 



(13) 



'rad.O 



10.33 
1.68 



(9-1) 



(14) 



where Hq is the present value of the Hubble constant and 
h is the Hubble constant in units of 100 km s _1 Mpc -1 . 
Without using the full solution for the matter-radiation 
era, existing here because of the validity of General Rel- 
ativity, it is possible to estimate the Hubble radius at 
equality by assuming that the matter solution holds all 
the way backwards in time. At equality, 

a P „H, 



a H 



= V2y/1 + z eq = 219/i [1 + (1 - <?)3.08] , (15) 
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where it was used a matter dominated solution for the 
scale factor in a flat Friedmann-Robertson- Walker uni- 
verse, a(t) oc i 2 / 3 , equivalently, the dependence of p m 
given m (12). The y/2 factor appears because at equal- 
ity, the right hand side of equation (11) has two equal 
energy densities. Again, the first term in the bracket is 
a complete standard result while the second one is the 
correction due to the change in the statistical model. 



C. Galaxy surveys and parameter estimation 

Now the analysis proceeds in much the same way as 
was done in Ref. for Brans-Dicke gravitation, which 
can be seen for details. The shift in the matter-radiation 
equality will lead to a shift in the maximum of the power 
spectrum of horizontal nature. Such a shift was consid- 
ered by Tegmark Jl6| who introduced a phenomenolog- 
ical parameter rj and considered how accurately it can 
be measured by SDSS. Tegmark gives two expected es- 
timates on Ary/?7, (see panel 3 and 4 of Fig. 1 of Ref. 
fL6||). The first one corresponds to the case in which 
all parameters are considered fixed -or measured by 
other experiments- except r\. The second one, to the 
case in which all parameters are to be extracted from 
SDSS alone. For the scale currently going non-linear 
(k ~ 0.1/iMpc" 1 ), Ar]/r] ~ 0.02 and Arj/rj ~ 0.1 respec- 



tively 
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The accuracy of the measurements on Arj/rj 
depends on the scale considered. The scale taken here 
-and also in |2l|- avoid one to worry about non-linearity 
or bias effects. These estimations can even be improved 
considering data of other satellites, like MAP or Planck, 
as explained in pl| . 

In the normal General Relativity situation, the focus 
is pointed towards the determination of f2 m /i, where f2 m 
is the matter density parameter. Changing Sl m h yields 
a horizontal shift to the power spectrum. In our sit- 
uation, using Arj/rj ~ 0.02 and equation (15) we can 
read limits upon q to be of order 10 -3 which is compara- 
ble to the nucleosynthesis bound obtained in [ll]|lj] and 
slightly less restrictive than the more in depth analysis 
of H and the bounds given in ju| and |H| . However, it 
comes now from a different arena, being the properties 
of galaxy surveys and the neutrino-photon temperature 
relationship its main inputs. 



D. Detectability of the cosmic neutrino background 

In the standard Boltzmann-Gibbs statistics, the ratio 
of the energy density of neutrinos to that of photons is, 
as we have noted before, 



!>) = 7 -(f) 4/3 N v = 0.681. 
pJbg 8 W 



(16) 



e + — e~ annihilation is not entirely correct ^,£3|. If 
the neutrinos share in the heating somewhat, their en- 
ergy density is larger than the above quoted value. This 
can be modeled as if the number of neutrinos families N v 
is effectively bigger than 3. Together with QED correc- 
tions to the energy density of e ± and 7's |2J] , the increase 
in the effective number N v is, 



5N V = 0.04 to 0.05. 



(17) 



This leads to an increase of ~ 1% in the neutrino energy 
density. 

In the non-extensive context, we may think of formula 
(9) as an effective variation in 7V„, given by 



SN V 



15 8, 



7 4 7 - 34 - 7 °(9- 



1) 



6.10(9- 1). 



(18) 



However, recent studies have pointed out that the as- 
sumption that neutrinos decoupled completely before 



Then, if a direct detection of the cosmic neutrino back- 
ground would be avalaible, a straightforward comparison 
with ( [l8| ) may also yield a bound upon (q — 1). 

The precision detection of the neutrino energy density 
was considered by Lopez et al. |2f|. They analized the 
possible consequences that a small increase in p v would 
have in the cosmic microwave background. New mea- 
surements planned for the forthcoming satellites MAP 
and Planck lead them to show that the sensitivity will 
be so great that the neutrino energy density will be de- 
tectable by itself. As in the work by Tegmark, estimates 
come in two different ways. If all other parameters (for 
instance, baryon density f2s, Hubble constant H , slope 
of the primordial perturbations n, etc) are held fixed (as 
given by other experiments), the expected sensitivity in 
N v is given by Fig. 3 of Ref. p5|] . Going up to multipole 
moment I ~ 1000 and using only temperature anisotropy 
data, the expected sensitivity will be ~ 0.01 |25|. If the 
standard result N v = 3 is confirmed, one can then obtain 
-using (p"8|)- a bound upon (q — 1) of order xlO -3 . This 
could even be improved using polarization data in the 
analysis, as explained in Ref. fl25j. 



E. Concluding remarks 

We have presented the case for obtaining new cos- 
mological bounds upon the degree of non-extensivity of 
the universe. These new bounds will come from accu- 
rate comological measurements of the microwave back- 
ground and galaxy surveys. To allow for this bounds to 
be obtained, we have computed the shift in the epoch 
of matter-radiation equality due to a different statistical 
background. Afterwards, comparing with the expected 
measurements of Arj/rj, we have shown how this can be 
used to constrain the possible non-extensive parameter. 
We have also analized the implications that a precision 
detection of the cosmic neutrino background may have 
concerning the fixing of the statistical description and 
we have shown that another bound may arise from there. 
The order of magnitude of these two bounds is expected 
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to be \q— 1| < 10 -3 and can be compared with other pre- 
viously obtained cosmological constraints ||ll]Jl3|] . How- 
ever it is important to stress here that many phenomena 
can cause the same effect. For instance, in the galaxy 
survey oriented scheme, a change in Q m or h, a different 
number of masless species or a change in the theory of 
gravity to a scalar-tensor one have all the same output 
on a eq H eq . In the case of the correction to the effective 
number of neutrino species N v , it may well be the case 
that the 1% correction arising from the extra heating 
and QED effects be hidden by corrections arising from 
a value of q different from 1. To determine this, the 
former corrections need to be computed directly in the 
non-extensive framework, without assuming Boltzmann- 
Gibbs statistics as an initial hypothesis. It is disturb- 
ing to note that whereas there is only one way of being 
standard there are many different and arbitrary ways to 
deviate from that situation. Even more disturbing is to 
sec that different deviations may produce the same obser- 
vational output. Despite of these caveats, observational 
facts like the one discussed here, as well as those concern- 
ing nucleosynthesis, are in the way of determining if the 
usual statistics is the right description of the universe. 
Although it of course seems a rather safe assumption, 
and in particular values of |g — 1| > 10 -3 are completely 
discarded, there is certainly the need for more work to 
get a definite answer. 
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